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Abstract 



We examine the composition of the L°° norm with weakly convergent sequences of gradient 
fields associated with the homogenization of second order divergence form partial differential 
equations with measurable coefficients. Here the sequences of coefficients are chosen to model 
heterogeneous media and are piecewise constant and highly oscillatory. We identify local rep- 
resentation formulas that in the fine phase limit provide upper bounds on the limit superior of 
the i°° norms of gradient fields. The local representation formulas are expressed in terms of the 
weak limit of the gradient fields and local corrector problems. The upper bounds may diverge 
according to the presence of rough interfaces. We also consider the fine phase limits for layered 
microstructures and for sufficiently smooth periodic microsturctures. For these cases we are 
able to provide explicit local formulas for the limit of the L°° norms of the associated sequence 
of gradient fields. Local representation formulas for lower bounds are obtained for fields cor- 
responding to continuously graded periodic microstructures as well as for general sequences of 
oscillatory coefficients. The representation formulas are applied to problems of optimal material 
design. 

Keywords: L°° norms, Nonlinear composition, Weak limits, Material design, Homogenization. 
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1 Introduction 

Understanding the composition of nonlinear functionals with weakly convergent sequences is a central 
issue in the direct methods of the calculus of variations, homogenization theory and nonlinear partial 
differential equations. In this paper we discuss a composition motivated by problems of optimal 
design. To fix ideas consider a domain SI C Mr, d = 2, 3, partitioned into two measurable subsets ui 
and Sl/w. Define the piecewise constant coefficient of thermal conductivity taking the values al for 
iew and f3I for x £ Q/w by A(uj) — (axw + /3(1 — Xu))I- Here is the characteristic function of 
oj with Xu = 1 for points in oj and zero otherwise and / is the d x d identity matrix. Next consider a 
sequence of sets {oJ n }^ = i with indicator functions % Wn and the iJ 1 (SI) solutions u n of the boundary 
value problems u n = g on 9S1 with g £ H 1 ^ 2 (dfl) and 



for / g _ff~ 1 (Sl). The theory of homogenization [5], [35], [32] asserts that there is a subsequence of 
sets, not relabeled, and a matrix valued coefficient A H (x) £ L°°(Sl,M. dxd ) for which the sequence u n 
converges weakly in i? 1 (Sl) to u H g if 1 (SI) with u H — g for x £ 9S1 and 



div (A(w„)V« n ) = / 



(1.1) 



div (A H Wu H ) = f. 



(1.2) 
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The compositions of interest are given by the L°° norm taken over open subsets S C f2 and are of 
the form 

l|Vu„|| L oo (s) = esssup xe5 |Vu n (x)|, (1.3) 

||X-„Vu„||l-(s) and ||(1 -Xa,„)Vu„|| L =c (s) , (1.4) 
and we seek to understand the behavior of limits of the kind given by 

liminf ||x w »Vit n |U«>(s) and Iimsup ||x M „Vtt n |U°°(s). (1-5) 

n ^oo n — 

In this paper we provide examples and identify conditions for which it is possible to represent the 
limits of these compositions by local formulas expressed in terms of the weak limit Vu H . The 
representation formulas provide a multi-scale description useful for studying the composition. 

To illustrate the ideas we display local formulas in the context of periodic homogenization. The 
unit period cell for the microstructure is Y and we partition it into two sets P and Y/P. To fix 
ideas we assume the set P represents a single smooth particle, e.g. an ellipsoid. The union of all 
particles taken over all periods is denoted by u. The coefficient A(u>) is a periodic simple function 
defined on K d taking the value al in cj and j3I in M. d /u. On rescaling by 1/n, n — 1,2,... the 
set given by the union of rescaled particles taken over all rescaled periods is denoted by uj n and 
Xu n i x ) — Xu(nx). We consider the sequence of coefficients A{oj n ) restricted to f2 and the theory of 
periodic homogenization [3] , |37j delivers a constant matrix A H of effective properties given by the 
formula 

A% = J^A lk (y)P kj {y)dy (1.6) 
where P k j = d Xk (f>>{y) + S k j and <j>> are F-periodic Hl oc {R d ) solutions of the unit cell problems 

div(A(j/)(V<^(y) + e j )) = in Y, (1.7) 

where we have written A(y) — A(uS) — (ctXu(y) + (1 — Xu{y))P)I f° r U . It is well known that 
the associated energies taken over sets 5CC converge [38], [32], i.e., 

lim / A n Vu n ■ Vu n dx = / A H Vu H ■ Vu H dx 
n ^°°Js Js 

A{y)P(y)Vu H (x) ■ V/(i) dydx. (1.8) 

SxY 

In this paper we show that the analogous formulas hold for L°° norms and are given by the local 
representation formulas 

lim \\Xu n { x )V u n\\L<*>(S) = \\Xu{y)P(y)Vu H (x)\\ L ^ {SxY ), and (1.9) 
lim ||(1- Xa^aOJVunllL-fs) = ||(l-Xa,(tf))i , (l/)Vu H (a:)|| li » ( sxy5, (1-10) 
these formulas follow from Theorem 1421 



For general situations the question of finding local formulas is delicate as the solutions of (1.1) 
with measurable coefficients are nominally in H 1 ^) with gradients in L 2 (il,R d ). For sufficiently 
regular /, g, and fl, and in the absence of any other hypothesis on the coefficients, the theorems 
of Bojarski [4 , for problems in M 2 , and Meyers [3D], for problems in M. d , d > 2, guarantee that 
gradients belong to L p (fl,M. d ) for 2 < p < p' with p' depending on the aspect ratio j3/a. For the 
general case one can not expect p to be too large. The recent work of Faraco [T^] shows that for d = 2 
and for j3 = K > 1 and a = 1/K that there exist coefficients associated with sequences of layered 
configurations u) n made up of hierarchal laminations for which the sequence of gradients is bounded 
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in Lf oc (Q,R rf ) for p < p* = 2K/(K - 1) and is divergent in Lf oc (n,R d ) for p > p*. This precise 
value for was proposed earlier for sequences of laminated structures using physical arguments in 
the work of Milton [3T]. For measurable matrix valued coefficients A(x) € M 2x2 with eigenvalues in 
the interval [1/K,K] the same critical exponent p* = 2K/(K — 1) holds, this result also motivated 
by [3T] is shown earlier in the work of Lionetti and Nesi [5T] . 

With these general results in mind we display, in section 2, a set of upper bounds on the limit 
superior of the compositions (1.4) that hold with a minimal set of hypothesis on the sequence 
{ w n}S° = i- Here we assume only that the sets Lo n are Lebesgue measurable thus the upper bound may 
diverge to oo for cases when these sets have corners or cusps. The upper bound is given by a local 
representation formula expressed in terms of the weak limit Wu H . It is given by the limit superior 
of a sequence of L°° norms of local corrector problems driven by Vu H . For periodic microstructures 
the local correctors reduce to the well known solutions of the periodic cell problems associated with 
periodic homogenization [3J, [37]. In section 3 we provide a general set of sufficient conditions for 



which the limits (1.5) agree and are given by a local representation formula see, Theorem 3.2 As 



before this formula is given in terms of a limit of a sequence of L°° norms for solutions of local 
corrector problems driven by Vu H . From a physical perspective the local formula measures the 
amplification or diminution of the gradient Wu H by the local microstructure. Formulas of this type 
have been developed earlier in the context of upper and lower bounds for the linear case [25] . [26] . 
[2"T] and lower bounds for the nonlinear case [T5] . 

On the other hand when the boundary of the sets u> are sufficiently regular one easily constructs 
examples of coefficients A(cj) for which the gradients belong to L°°(n,M. d ). More systematic treat- 
ments developed in the work of Bonnitier and Vogelius [5] , Li and Vogelius [19] , and Li and Nirenberg 
[2Ti] describe generic classes of coefficients A(uS) for which gradients of solutions belong to L^ c (fl, M. d ). 
The earlier work of Chipot, Kinderlehrer and Vergara-Caffarelli [5] establish higher regularity for 



coefficients A(uj) associated with laminated configurations. In section 4.1 we apply the uniform 



convergence for simple laminates discovered in [5] to show that the sufficient conditions given by 
Theorem 3.2 hold. We apply this observation to obtain an explicit local formula for the limits of 
compositions of the L°° norm with weakly convergent sequences of gradient fields associated with 
layered microstructures. While in section [4~2] we use the higher regularity theory for smooth periodic 
microstructures developed in |20j to recover an explicit representation formula for the upper bound 
on the limit superior of compositions of the L°° norm with weakly convergent sequences of gradient 
fields associated with periodic microstructures. Lower bounds on the limit inferior are developed in 
section 5 that agree with the upper bounds and we recover explicit local formulas for the limits of 
compositions of the L°° norm with weakly convergent sequences of gradient fields associated with 
periodic microstructures. 

The L°° norm of the field gradient inside each component material (1.4) is of interest in ap- 
plications where it is used to describe the strength of a composite structure. Here the strength of 
a component material is described by a threshold value of the L°° norm of the gradient. If the 
L°° norm exceeds the threshold inside ui n then failure is initiated in that material and nonlinear 
phenomena such as plasticity and material degradation occur [IB], [31] • The design of composite 
structures to forestall eventual failure initiation is of central interest for aerospace applications [16] . 
For a given set of structural loads one seeks configurations uj that keep the local gradient field below 
the failure threshold inside each component material over as much of the structure as possible. As 
is usual in design problems of this sort the problem is most often ill-posed (see, e.g. [28]) and there 
is no best configuration uj. Instead one looks to identify sequences of configurations {uJn}^^ from 
which a nearly optimal configuration can be chosen. 

The work of Duysinx and Bendsoe [10] presents an insightful engineering approach to the problem 
of optimal design subject to constraints on the sup norm of the local stress inside a laminated 
material. The subsequent work of Lipton and Stuebner [22], [23], [24] develops the mathematical 
theory and provides numerical schemes for the design of continuously graded multi-phase clastic 
composites with constraints on the L°° norm of the local stress or strain inside each material. More 
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recent work by Carlos-Bellido, Donoso and Pedregal [TTJ provides the mathematical relaxation of 
the L°° gradient constrained design problem for two-phase heat conducting materials. The feature 
common to all of these problems is that they involve weakly convergent sequences of gradients and 
their composition with L°° norms of the type given by ( 1.3 1 and ( 1.4 1. Motivated by the applications 
we develop an explicit local representation formula for the lower bound on (1.5) for continuously 
graded periodic microstructures introduced for optimal design problems in |25j . [22] , |23j , see section 
5. A similar set of lower bounds have appeared earlier within the context of two-scale homogenization 
|27j . In section 6 we conclude the paper by outlining the connection between optimal design problems 
with L°° gradient constraints, local representation formulas, and the composition of the L°° norm 
with sequences of gradients. Last it is pointed out that the results presented here can be extended 
without modification to the system of linear elasticity. 



2 Mathematical background and upper bounds given by lo- 
cal representation formulas 

In this section we present upper bounds on the limit superior of sequences of L°° norms of gradient 
fields associated with G-convergent sequences of coefficient matrices. In what follows the coefficient 
matrices given by simple functions A(x) taking the finite set of values A\,A2, ■ ■ ■ ,An in the space 
of d x d positive definite symmetric matrices. Here no assumption on the sets w, where A(x) = Aj 
are made other than that they are Lebesgue measurable subsets of f2. 

We consider a sequence of coefficient matrices A n (x) = J2iLiXnAi- Here A n (x) = A% on the 
sets iu n and the corresponding indicator function \n takes the value x n = 1 on w n an d zero outside 
for i = 1,2, ... ,N, with YliLiXn = 1 on O. We suppose that the sequence {A n (x)} n ° =1 is G- 
convergent with a G-limit given by the positive definite dx d coefficient matrix A H (x). The G- limit 
is often referred to as the homogenized coefficient matrix. For completeness we recall the definition 
of G-convergence as presented in [32] : 

Definition 2.1. The sequence of matrices {A n (x)}^' =1 is said to G-converge to A (x) iff for every 
lu C fl with closure also contained in fl and for every f G the solutions ip n G Hq(uj) of 

- div {A n V<p n ) = f (2.1) 

converge weakly in Hq(lj) to the Hq(lu) solution ip H of 

- div (A H Vtp H ) = f. (2.2) 

G-convergence is a form of convergence for solution operators and its relation to other notions 
of operator convergence are provided in [55]. From a physical perspective each choice of right 



hand side / in (2.1) can be thought of as an experiment with the physical response given by the 
solution tf n of (2.1). The physical response of heterogeneous materials with coefficients belonging 
to a G-convergent sequence converge in Hq(oj) to that of the G-limit for every choice of sub-domain 
uj. For sequences of oscillatory periodic and strictly stationary, ergodic random coefficients the G- 
convergence is described by the more well known notions of homogenization theory [5J, [17] . [36] . 



[55] . [52] . We point out that the G-convergence described in Definition 2.1 is a specialization of the 
notion of H-convergence introduced in [32] which applies to sequences of non-symmetric coefficient 
matrices subject to suitable coercivity and boundedness conditions. 

It is known [32; that if {A n }^ =1 G-converges to A H , then for any g G # 1/2 (<90) and / G H^iQ), 
the 7? 1 (fi) solutions u n of 

- div {A n Vu n ) = /, u n = g (2.3) 
converge weakly in (SI) to the iJ 1 (£7) solution u H of 

- div (A H Vu H ) = /, u H = g. (2.4) 
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Last we recall the sequential compactness property of G-convergence [38] , [32] applied to the case 
at hand. 

Theorem 2.1. Given any sequence of simple matrix valued junctions {A n (x)}^ ) =1 there exists a 
subsequence {A n (x)}^ =1 and a positive definite d x d matrix valued function A H (x) such that the 
sequence {A n (x)}^) =1 G-converges to A H (x). 

For the remainder of the paper we will suppose that sequence of coefficients {A n }^L 1 G-converges 
to A H and we will investigate the behavior of the gradient fields inside each of the sets uj l n . To this 
end we will consider the limits 




where S C f2 is an open set of interest with closure contained inside f2. 

In order to proceed we introduce the local corrector functions associated with the sequence of 
coefficients {A n }^ =1 . Let Y C K d be the unit cube centered at the origin. For r > consider 
fi* n * = {x £ fi : dist(x,dtt) > r} and for x £ f2* nt and z £ Y we introduce the Y periodic H 1 ^) 
solution w r,n (x, z) of 

- div z (A n (x + rz)(V z w^ n {x, z) + e)) = 0, for z £ Y, (2.6) 

where E is a constant vector in K d with respect to the z variable. Here x appears as a parameter 
and the differential operators with respect to the z variable are indicated by subscripts. For future 
reference we note that w r,n depends linearly on e and we define the corrector matrix P r,n (x,z) to 
be given by 

P r ' n (x,z)e = V z w^ n (x,z) + e. (2.7) 

We are interested in the L°° norm associated with each phase and introduce the modulation functions 
M l {Vu H ) defined for x £ fl given by [2"5] 

M i (Vu H )(x) =\misup\imsup\\xi(x + rz)(P r > n (x,z)Vu H {x))\\ L o. {Y) . (2.8) 
r— >0 n— >oo 

In what follows we will denote the measure of uj C f2 by and state the following upper bound 
given by a local representation formula. 

Theorem 2.2. Let A n G-converge to A H and consider any open set S C f2 with closure contained 
inside Q. There exists a subsequence, not relabeled and a sequence of decreasing measurable sets 
E n C S, with \E n \ \ such that 

\imsvL V \\xiyu n \\ L ~ iS \ Ert) <\\M l (Vu H )\\ L ~ {S ), i = l,2,...,N. (2.9) 

To proceed we introduce the distribution functions associated with the following sets S™ t , i = 
1,2, ...,N, defined by 

Sl t = {x £ S : xi\Vu n \ > t] (2.10) 

given by 

K^) = \Sl t \. (2.11) 
We state a second upper bound that follows from the homogenization constraint [25 . 
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Theorem 2.3. Let A n G-converge to A H and consider any open set S C O with closure contained 
inside Q. Suppose for i = 1, 2, . . . ,N that limsup^^^ ||x* l Vw„||ioo(5) = i l < oo and for every 6 > 
sufficiently small there exist positive numbers 9\ > such that 



liminfA?(f -8) > 0j. 

n— >oo 

There exists a subsequence, not relabeled, such that 



(2.12) 



(2.13) 



We provide a proof Theorem 2.2 noting that the proof of Theorem 2.3 



is given m 



Proof. First note that the claim holds trivially if ||.M l (Vu )||£°°(,g) = oo. Now suppose otherwise 
and set ||A^ I (VM ff )||ioo(5) = H < oo. For this case Corollary 3.3 of [55] shows directly that for any 
6 > that the measure of the sets 



S 



{xeS: xi(x)\Vu n (x)\>H + S}, 



tends to zero as n goes to oo, i.e., 

limsup \1(H + S)= limsup \S? H+S \ = 0. 



(2.14) 



(2.15) 



We choose a sequence of decreasing positive numbers {8i}fL ll such that Si \ and from (2.15) we 
can pick a subsequence of coefficients {A n ^ Sl ^ , }'^L l for which 
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71j(<5l) 

i,H+$i 



<2- J , j = 1,2,. 



(2.16) 



For 82 we appeal again to (2.15 1 and pick out a subsequence of {A n ^^ Sl '}^. 1 denoted by {A 7lj ^}°^ 1 



for which 
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n j (S 2 ) 
i,H+S 2 



<2-\ J = 1,2,.... 



(2.17) 



We repeat this process for each 81 to obtain a family of subsequences {A n i( s ^} S jL 1 , I = 1, 2, . . . such 
that {A n i {s *+^}f =l C {A n ^ s ^}f =l . On choosing the diagonal sequence {A nki - S ")}^ =1 we form the 
sets 



E K = U k > K Slf + % = {x&S: X \ 



k (8 k )\^ U n k (S k ) 



> H + 8 k , for some k > K}, (2.18) 



with E K+ i C E K . Noting that |5"^+/J < 2_fc > we see that \ e k\ < 2~ K+1 - Since x £ E K implies 
that 



Xn k (S k )\V U n k (5 k ) 



<H + 8 k for all k> K, 



we observe that 



\\Xn k (S k )^ u n k (5 k )\\L-(S\E K ) < H + 8 k for all k > K, 



and we conclude that 



limsup \\Xn K {SK) Vu nK{Su)\\L°°{S\E K ) < H, 
K— >oo 



with \Ek\ \ and the theorem is proved. 



(2.19) 
(2.20) 

(2.21) 

□ 
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3 Lower bounds and sufficient conditions for a local repre- 
sentation formula 

We suppose that sequence of coefficients {A n )}^ =1 G-converges to A H and investigate the behavior 
of the gradient fields inside each of the sets w^. Here we consider the limits 

liminf Hj&Vunlli-^, fori -1,2..., N. (3.1) 

and identify a general sufficient condition for obtaining a lower bound on these quantities in terms 
of M*{Vu B ). 

Assume that u n , P r,n and u H are defined as in the previous section and we consider an open 
subset S C ft with closure contained in f2. We write r = dist(dS, dfl) > and set 

S T = {x e £1 : dist(x, S) < r}. 

For r < t note that S C S r C 5 r C fi. We next recall for x £ S 

M\Wu H )(x) = \imsnp\imsup\\\xl % (x + ry)P^\x,y)Wu H (x)\ 2 \\ La , {Y} . (3.2) 

7 >-0 71— ^OO 

For this case the sufficient condition is based on the distribution function for the sequence {Xn(x+ 
ry)P r ' n (x,y)\7u H (x)} and the lower bound is presented in the following theorem. 

Theorem 3.1. Let A n G-converge to A H and consider any open set S C O with closure contained 
inside fi. Suppose that 

\\M l {Vu H )\\ L ~ [s) =t <oo. 
Assume also that for all S > small, there exist /3s > swc/i that 

lim liminf \{(x,y) e S X Y : |xU* + n/)P r ' n (z,y)Vu H (z)| 2 > (^) 2 - <J}| > fa > 0. (3.3) 

r— >0 n— ^oo 

T/ien i/iere exists a subsequence for which 

lim liminf ||x^VM„|| L = (Sr) > WM^Vu 11 )^^^. 

Proof. Our starting point is Lemma 5.5 of [7] which is described in the following lemma. 
Lemma 3.1. 

lim lim sup / / \P r ' n (x,y)Vu H (x) - V u n (x + ry)\ 2 dydx = 0. (3.4) 

On applying the lemma we observe that 
Xi(x + ry)P r > n (x, y)Vu H {x) = x\ l {x + ry)Wu n (x + ry) + z r > n (x, y) V(x, y) E S x Y, (3.5) 
where 

lim lim sup / / \z r > n {x, y)\ 2 dydx = 0. (3.6) 

n—>-oo Js JY 

It follows that 

\xUx + ry)P r ' n (x, y)\7u H (x)\ 2 = \ X l n {x + ry)Vu n (x + ry)\ 2 + F r ' n (x, y) V(x, y) e S x Y, (3.7) 
where 

F r ' n (x, y) = \z r > n (x, y)\ 2 + (z r ' n (x, y), tfjx + ry)Vu n {x + ry)). 
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We show that F r ' n (x,y) — > strongly in L 1 in the sense that 



lim lim sup 



r->0 



\F r ' n (x,y)\dydx = 0. 



n— too J 5 Jy 

Indeed, from the definition and by Cauchy-Schwarz inequality we have 



\F r > n (x,y)\dydx < 



S JY 



\z r > n (x,y)\ 2 dydx 



jJ Y \z r - n (x,y)\ 2 dydxy (JJ Y \Vu n 



(x + ry))\ 2 dydx 



1/2 



(3.8) 



Moreover from standard a-priori estimates we know there is a constant C > independent of r and 
n for which, 



\Vu n {x + ry))\ 2 dydx < C. 



(3.9) 



S JY 



The assertion follows on taking the limits in (3.8 1 and using estimate (3.9) and equation (3.6). 
Now by Chebyshev's inequality, for every S > 0, we have the inequality 

1 



\{(x,y)eSxY:\F r ' n (x,y)\>5}\< 



SxY 



\F r - n (x,y)\dydx 



(3.10) 



and taking the limsup as n — > oo first and then as r — > 0, we see that 

lim lim sup | {(x,y) G S x Y : \F r > n (x,y)\ > S}\ = 

From ( |3.7[ ) we see that 

{(x,y) eSxY: \^ n {x + ry)P r - n (x, y)Vu H {x)\ 2 > (tf - 6} C 

C {(x, y)eSxY: \ X l n (x + ry)Vu n {x + ry)\ 2 > (t) 2 - 26} U {(x, y)eSxY: \F r ' n (x, y)\ > 6}. 



Therefore, applying (3.10) we obtain 

hm liminf \{(x,y) eSxY : \ X Ux + ry)P r > n {x, y)Vu H {x)\ 2 > {tf - S}\ < 

r— >0 n—>oo 

< lim liminf \{(x, y) G S x Y : \yt{x + ry)Vu n {x + ry)\ 2 > (t f - 26}\. 

r— >-0 n— >co 

It follows from the last inequality that 

Km liminf | {(x,y) eSxY : \xi(x + ry)Vu n {x + ry)\ 2 > {tf - 25}\ > p s > 0, 

r— >0 n— »oo 



Here we have used our assumption ( 3.3). Therefore, there exist R = R{6) and N = N(6) such that 
\{(x,y) eSxY : \j&(x + ry)Wu n (x + ry)\ 2 > (t) 2 - 26}\ > 0, Vn > N(S),r < R(S). 



From the definition of the L°° norm it follows that 



Hlx^Vu'YlU^) > (ty - 26 Vn > N(5),r < R(6). 
Taking the limit first in n and then in r, and using the arbitrariness of of <5, we get 



lim liminf \\xl^u n \\l^ ( s r) > \ 



r— >0 n— >-oo 



and the theorem follows. 



□ 
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Last if we combine the hypotheses of theorems |2.3| and |3.1| we obtain a sufficient condition for 
a local representation formula for limits of compositions of the L°° norm with weakly convergent 
sequences of gradients associated with homogenization. 

Theorem 3.2. Let A n G-converge to A H and consider any open set S C O with closure contained 
inside Q. Suppose for sufficiently small r < t, S C S2r C S r C fl, for i = 1,2, . . . , N that 
lim sup,,^ o li m su Pn->oo l!Xn^ u n||L°°(s r ) = ( l < oo and for every d > sufficiently small there exist 
positive numbers 9\ > such that 

limsuplimsup|{a; G S r : xjjVu n | >t - 5)}\ > Q\ > 0, (3.11) 

r— >0 n— >oo 

in addition suppose that limsup r _ > , ||.M l (Vu )\\L°°(s r ) — £ l < oo and for all S > small, there exist 
Ps > such that 

limliminf \{{x,y) eS r xY: \^ n {x + ry)P r ' n (x, y)Vu H (x)| 2 > (f ) 2 - <5}| > ft > 0. (3.12) 

r— >0 n— >oo 

There exists a subsequence, not relabeled, such that 

lira lim ||x^V M „|| L oo (Sr) = lim \\M l (Vu H )\\ L ~ {Sr) . (3.13) 



Local representation formula for layered and periodic mi- 
crostructures 



We now describe sequences of configurations for which one has equality in the spirit of (3.13|. The 
first class of configurations are given by sequences of finely layered media. The second class is given 
by a sequence of progressively finer periodic microstructures comprised of inclusions with smooth 



boundaries. In what follows the results of [8] provide the sufficient conditions (3.11) and (3.12) for 
the case of finely layered media. While the higher regularity results of [19] and [20] allow for the 
computation of an upper bound for the periodic case. This upper bound agrees with an explicit 
lower bound developed in section 5. We note that the lower bound for the periodic case can also be 
obtained using the earlier results given in |27j . 

In order to proceed let us recall the fundamental results from homogenization theory for periodic 
media. We denote a d dimensional cube centered at x and of side length r by Y(x,r). For the 
unit cube centered at the origin we abbreviate the notation and write Y. The coefficient A(y) is 
a periodic simple function defined on the unit period cell Y taking the N values Ai, i = 1, . . . ,N 
in the space of positive symmetric d x d matrices. We denote the indicator functions of the sets 
Yi where A(y) = Ai by \ l an d write A(y) = A-iX l {y)- It is well known from the theory of 

periodic homogenization [3] that the sequence of coefficients A n (x) = A(nx) G— converge to the 
homogenized constant matrix A H given by the formula 



A% = j^A lk {y)P kj {y)dy (4.1) 

where P k j — d Xk <P{y) + 6 k j and <fp are F-periodic Hl oc (R d ) solutions of the cell problems 

div(A(y)(V^'(y) + e>)) = in R d , (4.2) 
where this equation is understood in the weak sense, i.e., 

(A(y) (V^(y) + e»'), Vii)dy - 0, W> G H l per {Y). (4.3) 

For periodic microstructures we define the modulation function by 

M\Vu H )(x) = \\\ X l (-)P(-)Vu H (x)\\\ L ^ {Y) i = l,..-,N. (4.4) 



Y 
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4.1 Laminated microstructure 



The layered configurations as introduced in this section are a special class of periodic configurations. 
To fix ideas we consider a two dimensional problem and partition the unit period square FcK 2 for 
the layered material as follows: 



Y 1 ={(y l ,y 2 )eY;-^<y 1 <-~ 



Y 2 = {(y u y 2 )eY:- 



1 



where 9 is a specified value in the interval (0, 1). Let x 1 and x 2 denote the indicator functions of Y\ 
and Y 2 respectively and consider the Y-periodic matrix function A(y) given by 

A{y) = a i x 1 (y) + (iix 2 {y), 

for positive constants a < (3. I is the 2x2 identity matrix. Let C K 2 and u n be the 
solution to 

- div (A(nx)Vu n ) = f in Q and u n = on dQ. (4.5) 
Then u n converges weakly in 1 (SI) as n — >• oo to the _ff 1 (il) solution u H of 



div (A H V 



H*-,„.H\ 



/, in O and u — on dfl 



(4.6) 



where A H is determined using the formula (4.1 1. The gradient of solutions of the cell problem (4.2) 
for layered materials are given by 



V*'(|f) 



and 



We define the constants 
ah = 



1 + (l-6)a 
Vcj) 2 {y) = e 2 for ally eY 

a/3 



+ (l-#)a X W J 



and a r . 



0a + (1-0)0, 



(4.7) 



0,5 + (1 - 0)a 

and introduce the Y periodic scalar coefficient a(y) — a\ 1 {y) + Px 2 {v)- A simple calculation gives 

ft/, 



pu(y) o 
o 1 



where pii(y) 



o(y) 



The homogenized matrix A ff is given by 



.4 



H 



a h 
a m 



The modulation function for each phase is given by 
M x {Vu H ){x) = 



-d u H 
6p+{l-6)a Xl 



M 2 {Vu H )(x) 



0/3+(l-0)a 



d Xl u*\ +(d X2 u H y 



We now apply the regularity and convergence results associated with G-convergent coefficients for 
sequences of layered materials |S]. For right hand sides / € if 1 (0) there exists a p > 2 such that for 
any subdomain Q' <e Q 

u n eil 1 ' 00 ^') and d X2 u n , a{nx)d Xl u n £ H^ P (Q!) 
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with the estimate that for some C = C(a, j3, O', f2), 

\\dx 2 u n \\H^p(n>) + \\a{nx)d Xl u n \\ H i, V {n>) < C\\f\\ s i(a), 



(4.8) 



see [8]. The Sobolev embedding theorem implies that {5 a:2 u„}^ 1 and {a(nx)d Xl u n }^ =1 are equicon- 
tinuous families over f2' and uniformly bounded in C(fi'). Then from ( |4.8| and the weak convergence 
U„ — * u H in H l {£l) it follows that for a subsequence 

d X2 u n — > d X2 u H , a(nx)d Xl u n — > cihd Xl u H uniformly in fl'. (4-9) 

We observe that 

al^iin -pn(na;)9 Xl w H | < a(nx)|a 2;i u„ - p 11 (nx)d Xl u H | = ^(na;)^^ -a^t^u^l, (4.10) 



and on applying (4.9) and noting that P(y) is constant inside each phase we see for i = 1,2 that 

|x l (nx)Vu n | = |x l (nx)P(na;)Vw ff | +m;(a;) (4.11) 
= >P(Vti H )(a:) + r<(x) (4.12) 

where m^(x) — » uniformly in f2'. Hence we arrive at the local representation formula for layered 
microstuctures given by 

Theorem 4.1. 

hm \\x l (nx)Wu n \\ Laa{ni) = ||MVu")|U-(no- (4.13) 

n— too 

It is easily seen that the uniform convergence implies that sequence of the gradients{Vu„} satisfy 
the non-concentrating conditions given by (3.11). Indeed, setting 

L l = limn-^co ||x I ( na -)^' u ™||i 00 (n')i an d for any 5 > there exists sufficiently large n for which 
|mj,(a;)| < § for ieff and 

\ X l (nx)Vu n (x)\ > \M l Vu H (x)\- 5 - 

so 



{x en' : M\Vu H {x)) > L l - -} C {xen' :\ X i {nx)Vu n \> V -5}. 
Therefore we conclude that for L l > 5 > 

c 

liminf \{x G n' : |x<(nx)Vu n | > U — -}| > |{x G O' : X ,; (Vu ff )(a;) > L ,: - <S}| > 0. 

n— foo 2 



Last the non-concentrating condition (3.12) follows immediately from the piecewise constant nature 
of the corrector matrix P(y) for layered materials. 



4.2 Periodic microstructure 

We consider periodic microstructures associated with particle and fiber reinforced composites. As 
before we divide Y into a union of iV disjoint subdomains Y\ . . . Y/v. Instead of proceeding within the 
general context developed in [20], [19] we fix ideas we suppose that the domains Y±, . . . , Y}v-i denote 
convex particles with smooth (i.e., C 2 ) boundaries embedded inside a connected phase described by 
the domain Y/v, se e Figure [TJ As before we denote the indicator function of Yj by x % an d the Y 
periodic coefficient is written A(y) = X l {v)Ai with each Ai being a symmetric d x d matrix of 

constants satisfying the coercivity and boundedness conditions given by 

A|£| 2 < (M, < A|C| 2 G K rf , and i = 1, . . . , N. 
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Figure 1: Particle reinforced geometry for two inclusions Y\ and Y 2 . 



For any bounded domain £1 C R d we consider the H 1 ^) solutions u n of 

div(A(nx)S/u n ) =0 in fl 



(4.14) 



associated with prescribed Neumann or Dirichlet boundary conditions. From the theory of periodic 
homogenization the solutions converge weakly in H 1 to the homogenized solution u . In this section 
we establish the following local representation theorem. 



Theorem 4.2. Let A(y) and the subdomains {Y i }fL 1 be as described above. Suppose u n solves (4.14) 
and u H is the corresponding homogenized solution, then for any subdomain fi' compactly contained 
inside fl one has the local representation formula given by 



lim \\x l (nx)Vu n \\ L ^ {(V)) = \\M l (Vu H )\\ 



L°°(n')- 



(4.15) 



For the proof we will use the W 1,00 estimate for weak solutions of linear equations with oscillatory 
periodic coefficients obtained in [2] for smooth coefficients and later extended in [20] to include 
discontinuous but locally Holder coefficients. A W 1,p estimate for p < oo is given in [6j. We 
point out that we have restricted the discussion to periodic homogenization for particle reinforced 
configurations of the kind illustrated in Figure [T] However the regularity theory for oscillatory 
periodic coefficients developed in (20] applies to more general types of domains Yj , . . . , Y/v with 
C 1,a boundaries. We note that the proof given here goes through verbatim for period cells with 
coefficients satisfying the general hypotheses described in [2"U] . 

Theorem 1.9 of [2D] and a suitable rescaling shows that for any r > and Y(xq,t) C VI that 
there exists a positive constant C independent of r and n for which 

||VUn|U°°(y(a;o,r/2)) < Cr^ 1 \\u n \\ L ^ {Y{x ^ r)) . (4.16) 

The local L°° estimate for weak solutions of elliptic linear problems (Theorem 8.17, [M]) gives 

|Wll/»(y(x ,r)) < Cr ~ d/2 \\ u n\\L 2 (Y(x ,2r)), (4-17) 

where the constant C is independent of n and r. Combining the two estimates delivers the following 
lemma. 

Lemma 4.1. Let A(y) and the subdomains {Y{\ be as described above. Choose r € (0, 1) such that 



Y(xo,2r) C fi. Then if u n solves (4.14), then there exists C, independent of n and r such that 

\\Vu n \\ L =a(Y(x ,r/2)) < Cr^^ \\u n \\ L 2 (Y ( Xo ,2r)) ■ (4-18) 
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Proof of Theorem \4-S\ To prove the theorem we first show that there is a subsequence, rik, for which 

lim \\Vu nk (x) - P{n k x)\7u H {x)\\ L ~ i(V )- (4.19) 

k— foo 

To begin we choose xq € f2 and r > such that rn is an integer and Y(xq, r) d D, contains an integral 
number of periods of diameter 1/n. Then from (4.2) we see that (1 / 'n)4> 3 1 (n) is a Y(xq, r)-periodic 
H\ oc function satisfying 

div(A{nx)(\7(-(f>J(nx)) + e 1 )) = in R d (4.20) 
n 

Combining equations (4.14) and ( ]4.20 1 we note that 

div(A(nx)[Vu n - (Vw n {x,x Q ) + Vu H (x ))}) = in Y(x Ql r) (4.21) 

where 

d 1 

w n (x,x ) = V* -<fp (nx))d x .u H (x ) + u h (x ). 
z — < n 

j 

Observe that for this choice of Y(xo,r) 

Vw n (x, Xq) + X7u H (xo) = P(nx)Vu H (xq) 

Adding and subtracting P(nx)Vu H (xq) delivers 

||Vu„(a;) - P{nx)\Ju H (x)\\ L ^(Y( XQ ^ r /2))) < ||Vu n (x) - P(nx)Vu H {x )\\ l ^(y( Xq ,t/2))) 

+ \\P(nx)Vu H (x) - P(nx)Vu H (x )\\ L 2 {Y{xo , r/2))) . (4.22) 



We apply Lemma |4.1| to find a constant C independent of n and r such that the following estimate 
holds true: 

\\Vu n -P(nx)Vu H (xo)\\ L °°(Y(x ,r/2))) 
C 

;\\u n - (w„ + V/(x ) • (x - x ))\\ L 2 (Y(x(ur))) 



< 



r -(d+2)/2 I 



Combining with (4.22) we obtain 



\\Vu n -P{nx)Vu H {x)\\ L ~> i Y(x ,r/2))) 

c 

r - (d+2)/2 H M » ~ (w n + Vu H (x ) ■ (x-X ))\\ L 2(Y( Xo , r ))) 



< 



(4.23) 



+ \\P(nx)Vu H (x) - P(nx)Vu H (xo)\\ L 2(Y(x ,r)))- 



We bound the the first and second terms on the righthand side of (4.23). The first term in the right 
hand side of (4.23) is bounded above by 

\\u n ~ (w n + Wu H (x Q ) ■ (x - X ))\\ L 2( Y ( Xo ,r))) 

< \\u H (x) - (u H (x Q ) + Wu H (x Q ) ■ (x - xo))\\ L 2 {Y(xo r)) 



1 

+ \\u n - U H \\ L 2(Y(x ,r)) + II ^2 - ft ( nx )dx 3 U H (x )\\ L 2 (Y 



(4.24) 



We apply Lemma |4.1| together with a priori elliptic estimates to find that 

||V^'(n-)|U~ (Kt n < C, 



(4.25) 
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where C is independent of n. Moreover, as u H is a solution of a PDE in divergence form with 
constant coefficients it satisfies 



\u H {x) - (u H (x Q ) + V/(r ) • (x - x ))\ < M\x - x \ 2 , x € ft' 
\Vu H {x) - Vu H (x )\ < M\x -x \, xeff 

where M is the supremum of \D 2 u H (x)\ over O'. Applying | ]4.25 1 and (4.26) gives 

\\Un - K + V!l H (l ) • (X - X ))\\ L 2 {Y ( Xo ,r))) 



(4.26) 



< C (r 2+d / 2 M + -L + \\u n - u H \\ L 2 {YM) 



(4.27) 



for some constant C independent of r and n. From periodicity it follows that 

\\V4P{nx)\\ L ^ [Y ( x ^r)) = \W^{v)\\l^{y) 



and applying Lemma 4.1 together with (4.261 delivers 

\\P{nx)Vu H (x) - P{nx)Vu H (x Q )\\ L 2 {Y{xQir)) < cA d+2 ^ 2 (4.28) 
From the theory of periodic homogenization see, [3J, |17) . one has the convergence rate given by 

\\u n - uH \\l 2 (Y(x ,t)) < C—j= 

V n 

and collecting results we have 

||Vu„- P(nx)Vu H {x)\\ L ~ {Y{x0tr/m < C (Mr + ^JL-) . (4.29) 



We pass to a subsequence rik, and consider Y(xo,rk) such that, r k — > 0, r k n k is an integer, and 
r[ d+2 — > oo as k — > oo. Then consider any subdomain fi' CC fl, and cover it with cubes 
{Y(xi, 7'fc/2)} a ; i6 o' . Using compactness we choose finitely many cubes so that 

fi'cu£ =1 y(a*,r fc /2), 

Now since ||Vu„ fc — P(n k x)Vu H (x)\\Loo( n t^ is bounded above by the L°° norms over a finite collection 
of cubes, we see that 



||Vn n , - P{n k x)Vu H {x)\\ Lx{ni)) <C\r k + — 1 



for sufficiently large k to conclude 

lim ||VM n J|£°°(n')) = lim \\P(n k x)Vu H (x)\\ L co^^, 



k—toc 



k—too 



(4.30) 



lim \\x l (n k x)Wu nk \\ L o, (n ,y } = lim^ \\x l (n k x)P(n k x)Wu H (x)\\ Lx(ni)) . 



(4.31) 



Now we bound (4.31) from above and below by ||.M*(V« )||i,°°(n')- First note for each n k and 
ie!!' that 



\ X l (n k x)P(n k x)Vu H {x)\ < \\ X l (-)P(-)Vu H (x)\\ L ~ (Y) 



(4.32) 
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and we conclude that 



lim \\x l (n k x)P(n k x)Vu H \x)\\ L «, m < ||Af (Vu H )|Uao (n ,). (4.33) 

k— >oo 



The lower bound 

||M*(V« H )|| z - ( no < lim ||x 4 Kx)V U „J| L oo (o0) (4.34) 
follows from a direct application of Corollary |5.1| proved the next section and we conclude that 

(4.35) 

The theorem follows on noting that identical arguments can be applied to every subsequence of 
{X l (nx)Vu n }^' =1 to conclude the existence of a further subsequence with limit given by (4.35). 

□ 



lim ||x l («-fc2:)Vu„JjL~(f2')) = \\M l (Vu H )\\ L °a( n ,). 

k— >OQ 



5 Continuously graded microstructures 

In this section we consider a class of coefficient matrices associated with continuously graded com- 
posites made from N distinct materials. In order to express the continuous gradation of the mi- 
crostructure we introduce the characteristic functions x z ( x j 2/) , i = 1 , . . . , iV belonging to L 1 (f2 x Y) 
such that for each x the function x l { x i is periodic and represents the characteristic function of the 
i th material inside the unit period cell Y. The characteristic functions are taken to be continuous 
in the x variable according to the following continuity condition given by 

lim / \x l (x + h,y)- x {x,y)\dy = 0. (5.1) 

The coefficient associated with each material is denoted by A4 and is a constant symmetric matrix 
satisfying the ellipticity condition 

A < A, < A 

for fixed positive numbers A < A. Wc define the coefficient matrix 

N 

A{x,y)= s }2A t x\x,y). 

i=l 

This type of coefficient matrix appears in prototypical problems where one seeks to design structural 
components made from functionally graded materials and [35]. Here the configuration of the 
N materials is locally periodic but changes across the domain Q. The composite is constructed by 
dividing the domain Q into subdomains flk.i, 1 = 1,... ,Mk of diameter less than or equal to 1/fc, 
k = 1,2,... and ft — uf^flk.i- Each subdomain contains a periodic configuration of N materials. 
The following lemma allows us to approximate the ideal continuously graded material by a piecewise 
periodic functionally graded material that can be manufactured. 

Lemma 5.1. For a given subdivision Qk,ii ■ ■ ■ &k,M k °f diameter less than 1/fc and any i = 1, • • • , N , 
there exists a sequence {Xk( x i v )}k^i of approximations to x l (x,y) given by 

I 

with the property that 

, lim / \xl( x >y)-x*(%,y)\dydx = o. (5.3) 
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In (5.2), xo fc i {x) denotes the characteristic function of £lk,i o,nd x\ i{y) — X l ( x k,i,y) is the char- 
acteristic function associated with the configuration of the i th phase inside the subdomain il^j °^ a 
sample point x^j S £lk,i- 



Proof. The definition of the approximating function is given in ( 5.2 1 . We verify that ( 5.3 ) is satisfied. 
For each x G fl, define the sequence of functions 



\x\{x,y) - x l {x,y)\dy. 



Then T l k (x) for all x € f2. Indeed, for a fixed x € f2, there exists a sequence of subdomains 
x € f2fc and points Xk,i k € &k,i k such that by definition, 



\x % (vk,l h ,y) - X l {x,y)\dy. 



It is evident that \xk,i h — x\ < 1/k since Xk,i k and x both belong to £)>k,l k - Applying the continuity 



condition (5.1), we see that ri(x) — >• as k 



convergence theorem. 



oo and (5.3) follows from the Lebesgue dominated 

□ 



Let us define the coefficient matrix of the functionally graded material. Divide the domain 
f2 into subdomains Clk,i, I — 1, ••• ) Affc of diameter less than or equal to 1/k, k = 1,2,... and 



£1 = (J l= \flk,i- Each subdomain contains a periodic configuration of N materials with period 1/n 



such that 1/k > 1/n. The configuration of the i th phase inside a functionally graded composite 
is described by x%( x i nx )> where x\( x i1J) ^ s gi ven by 
denoted by A k (x, nx) and is written as 



5.2 



The corresponding coefficient matrix is 



Mk 



A k (x, nx) = V xi(x, nx)Ai 



(5.4) 



As seen from the proof of the lemma the continuity condition ( 5.1 ) insures that near by subdomains 
£lk,i and £lk,l' have configurations that are nearly the same when 1/k is sufficiently small. The fine- 
scale limit of such composites is obtained by considering a family of partitions indexed by j = 1,2,..., 

k ■ 

with subdomains Q l 3 of diameter less that or equal to 1/ kj . The scale of the microstructure is given 
by 1 /n,j . Both 1 / kj and 1 /rij approach zero as j goes to infinity and we require that lim^oo jjk 
For future reference the associated indicator functions and coefficients are written 



0. 



Let 



Xk (x,njx) and A kj {x,njx). 



A H {x) = j^A(x,y)P(x,y)dy 



where the matrix P(x, y) is defined by 

n/ \ dw j . 
P{x,y) id = — + SH, 

and w l (x, •) is a Y periodic function that solves the PDE 

div y (A{x, y)(V v w\x, y) + e 1 )) = 0, 



(5.5) 
(5.6) 

(5.7) 
(5.8) 



where {e 1 }, i = 1, ... is an orthonormal basis for R d . 

The Sobolev space of square integrable functions with square integrable derivatives periodic on 



Y is denoted by H^Y). The functions w l (x,y) belong to C(fi; ^ er (Y)) this follows from JITT] ) 
and is proved in the Appendix. 

We present the homogenization theorem for the sequences A kj (x, njx) proved in 
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Lemma 5.2. One can construct sequences {x\ i x > xn j)}j°=i for which the coefficient matrices {A kj (x, njx)}j°^ 1 
G— converge to the effective tensor A (x) defined by (5.6). 

Let / £ H^{VL) be given. Then by Lemma 5.2 the sequence of solutions{wj } of the equation: 

-di\[A ki {x, rijx)Vuj{x)] = f, Uj £ Hq(CI) 

converge to u H weakly in Hq, where u solves the the equation 

~div[A E (x))Vu H (x)} = f, /g^fi). 

We now have the following result. 

Theorem 5.1. Let V C L l (Y) n L°°(Y) be a countable dense subset o/L 1 (F). Assume that all 
elements of V are periodically extended to M. d . Suppose that </>(x) £ C(f2), t](x) £ V and Uj, P and 
U be given as above. Then 

lim / (j)(x)r](n j x)xlAx,n j x)\Vu j \ 2 dx = / 4>(x)r](y)x l (x, y)\P{x, y)Vu H {x)\ 2 dydx. 

By taking the modulation functions for continuously graded composites (see, [25]) to be 

M\Vu H ){x) = \\x l {x, -)P(x, OVu^aOHj^y). (5.9) 
we obtain the following corollary. 
Corollary 5.1. 

||M*(Vu H )|| i/ oo( n) < limsup||xi i (or,rijx)Vu j ||i,«»(n). (5.10) 

j-yoo 

Under an additional asymptotic condition on the distribution functions for the sequence {Vuj}'^. 1 , 
equality can be achieved in the above corollary. Indeed, define 

Sli = {x:x%( x ,n j x)\Vu j (x)\ 2 >t}, x^(x) := x s i . (5.11) 

and the distribution functions are given by 

|S*,| = / xlidx. (5.12) 
Jn 

Passing to a subsequence, there exists density functions 9t^ such that 

Xi 4 (x) ± 6 t ,i(x) L°°weak* 

and for any open subset 5cfi 

lim \S 3 ti r\S\ = [ 9 tl dx. (5.13) 

We present a sufficient condition on the distribution functions |Sj J associated with {S7uj}°^ 1 for 
which equality holds in (5.10). 

Corollary 5.2. Suppose that I = lim sup^^ \\Xk- ( x > n j x )\^ u j\ ||i°°(n) < oo and for each S > 
there exists a positive number 0s > for which 

\{x£tt: 8^ > 0}| > S . (5.14) 

Then 

limsup||xfc j (x,n i x)Vu i || i c 0(n) = ||Af (V^H^o) 
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Proof. (Proof of Corollary 5.2) The homogenization constraint [25] states that for almost every 

x e n 

e t j(x)(M l {Vu H (x)) - t) > i = l,...N. 
It follows that on the set where 6 t ,i > 0, we have M. i (Vu H (x)) > t. Let 

h = llxy^^)|Vu-T||L~(n)- 
For a subsequence lj —> I. Then given S > 0, there exists a natural number J such that 
I - 5/2 < I, = || X 4 fcj (x,n,x)|W| 2 || L oo (0) < I + 6/2 Vj > J. 

The measure of the set S^_ s ^ 2 . is positive. Moreover Sf__ g ^ 2 ,. C Sf_ Si and 

Xls.i A 0l-S,i(x) L°° weak * as j -> oo 
From hypothesis the set where 9is,i > is a set of positive measure for all 6 > 0. Therefore, 

A<f(Vu H 0)) >i-<y, 

on a set of positive measure that is ||.<M I (Vif ff (x))|| OC) > / — 8. The corollary is proved since 8 > is 
arbitrary. □ 



Proof. (Proof of corollary 5.1 1 From Theorem 5.1 it follows that for any <f> g C(f2) and rj G V is V- 
periodic, 



X i (x,y)0(x)r ? (y)|P(a;,y)Vit ff | 2 d?/da; 



< lim / |^(x)r?(rijx)|£telimsup||xL(a:,Ti j a;)|Vuj| |U°°(fi) 
By the Riemann-Lebesgue lemma, 

lim / \<p(x)rj(n,jx)\dx = / 4>(x)dx / |^(i/)|o?j/. 
i^°oJn Jn Jy 

Dividing both sides by the L 1 -norm of 4>, we obtain that for every x G \ Z, where Z is a set of 
measure zero, 

X i (a;, 2 /)?7(?/)|P(x,y)Vw H | 2 d 2 / < / \r)(y)\dy lim sup ||xj. (x, rijx)| Vuj | 2 ||z~(fi) 

The set Z depends on the choice of rj. But since V is countable, the union of the sets Z corresponding 
to elements of V will be of measure zero and the above inequality is true for any rj £ V and for every 
x outside this union. Now divide the last inequality by the L 1 norm of r/ in Y. Taking the sup over 
V and noting that V is dense in proves the corollary. □ 



Proof. (Proof of Theorem 5.1 1 For j3 > , define 



Now let Vj solve 



Ax{x,y) = A{x,y) + Px l (x,y)<f>(x)v(y)I- (5-15) 
-div[A^ (x, rijx)Vvj] = /, Vj e Hq(Q) 
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Then for any <p € Hq{Q), we have 



(A^ (x, njx)Vvj, Vip)dx = / fifdx and 

£2 if! 

(A fej (x, n-jX)Vuj, Vip)dx — / /(^da; 
n Jn 

Let = Vj — i/j. Then subtracting the second equation above from the first, we obtain 

/ (A k 1 3 {x,n J x)V5u. ] ,Vtp)dx+ / {{A^ (x, rijx) - A k] (x, n i x))Vu :/ , V^)ofe = 0, 
in Jn 

for all € Jfo(f2). Simplifying the above equation we get 

/ {A\(x, fijx)V5uj, Vip)dx + (3 I x\- i x i Tijx)(f)(x)r](njx)(Vuj, Vip)dx — 0, 
Jn Jn 3 

Plug in tp = Uj in the above equation to get, 



(A^ (x, njx)W5uj, Wuj)dx + f3 / x\{ x i n j x )4>{ x )v{ n j x )W u j\ dx = 0, 
n Jn 3 



(5.16) 



(5.17) 



Also plugging in ip = Suj in (5.161 yields 



(A 3 (x, njx)Vuj, X75iij)dx = / fSujdx. 
n Jn 



(5.18) 



Subtracting (5.18) from (5.17) and noting that the coefficient matrices are symmetric we get 
Jn 3 Jn 

where 

T 3 = j3 I x\ 3 i x i njx)4>(x)r)(njx)(\7uj, VSuj)dx. 



Let us estimate T 3 . To begin with, observe that 

{A\ 3 (x, rijx)V 8uj, \75uj)dx + j3 \ x\ ( x ,n'jx)4>(x)Tj(njx)(yuj, VSiij)dx = 0, 
n Jn 3 

Them from ellipticity, we get 

a / | VSuj \ 2 dx < / (A^ 3 (x, njx)V5iij, VSuj)dx 
Jn Jn 

<P Xk 1 ( x ' n j x )\H x ) r l{njx)\\'7uj\\V6v, j \dx 
Jn 

< C*/3||V(5u J '|| L 2||Vu J || L 2. 



That is 



\V6uA\ L 2<C/3, 



since the sequence Vm- 5 is bounded in L 2 . From this and the definition on T 3 we obtain 

\T j \ < C(3 2 
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From Lemma 5.2 we know that Uj — 1 u H , and Vj — 1 where and satisfy the following 
equations, respectively: For any <p € ifg(f2) 



JQ JQ 



A H (x)Wu H ,W(p)dx = / /pcfcc 

Q JQ 



(5.19) 



where ^4 ff (x), is the effective matrix given by (5.6 1 and 



A?(x) = J^A 1 {x,y)P 1 (x,y)dy 
where the matrix Pi(x,y) is defined by 

and w\(x, •) is a Y periodic function that solves the PDE 

div y (A 1 {x,y)(V y w\(x,y) + e 1 )) = 0, 

where {e 1 }, i = 1, ... is an orthonormal basis for R d . 

Writing 8u H = v H — u H and letting j — > oo, we obtain 



(5.20) 



(5.21) 



(5.22) 



^°°Jq 



= - t fSu H dx. 

JQ 



(5.23) 



One easily verifies that the variational formulations (5.19) can be written in terms of the two scale 
variational principles [33], [T] given by 



A 1 (x,y)(\7v H (x) + V y Vi(x,y)), V<p(x) + V y tp 1 {x,y)dydx = / fipdx 



n Jy 



Q JY 



A(x, y)(Vu H (x) + V v ui(x, y)), Vip(x) + V y (fi(x,y)dydx = / fipdx 



(5.24) 



where the solutions (u H , u\), (v H , vi), and trial fields (ip, ipi) belong to the space Hq (Q)xL 2 [fl; W££(Y)]. 
On writing 8u\ = v\ — u\, 5u H = v H — u , Ai(x,y) = A(x,y) + j3x 1 {x,y)<t>{x)rj(y)I , substitution 
into the first equation in (5.24) and applying the second equation in (5.24) gives 



A 1 (x,y)(ySu H (x) + V y 5u!(x,y)), V<p(x) + V y ipx(x,y)dydx 



Q JY 



+ f f (x l (x,y)^Hy)(^u H (x)+V y u 1 (x,y)),Vp(x)+V yipi (x,y)dydx 



Q JY 



0. 



Next we substitute (ip, <p{) — (8u H ,8u\) into the second equation of (5.24) to obtain the identity 
A(x, y)(Vu H (x) + V y Ui(x, y)), VSu H (x) + V y 5u\(x, y)dydx 



Q JY 



f(x)6u H (x)dx. 



(5.26) 
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On choosing (ip,(pi) — (u H ,ui) in (5.251 and applying (5.26) we obtain 



T + P / {x l (x,ymx)r 1 {y){Vu H (x) + V v u 1 (x,y)),Vip(x) + V y ip 1 (x,y))dydx 
Jn Jy 

fSu H dx, 



where 



T = P / / { X % {x 1 y)(j3{x)T 1 {y){^u H {x) + ^ y u 1 (x,y)) ) V5u H {x) + ^ y 8u 1 {x ) y))dydx. 
Jn Jy 



Next we set (92,^1) = (Su H ,Sui) in (5.25) and applying ellipticity delivers the estimate 

\\VSu H + V y Su 1 \\ L 2 (nxY) < CP, 

and we find that 

|T| < CP 2 . 



(5.27) 



(5.28) 



(5.29) 



(5.30) 



Since (5.23) and (5.27) have the same right hand sides we equate them and the theorem follows on 
identifying like powers of p. □ 

We conclude this section noting that lower bounds similar to those given here can be obtained 
in the context of two-scale convergent coefficient matrices [27] . 



6 Local representation formulas and gradient constrained de- 
sign for graded materials 

In view of applications it is important to identify graded material properties that deliver a desired 
level of structural performance while at the same time provide a hedge against failure initiation |16) . 
In many applications there is a separation of scales and the material configurations forming up the 
microstructure exist on length scales significantly smaller than the characteristic length scale of the 
loading. Under this hypothesis the structural properties are modeled using effective thermophysical 
properties that depend upon features of the underlying micro-geometry, see [13] , [29j . In this context 
overall structural performance measured by resonance frequency and structural stiffness are recovered 
from the solutions of homogenized equations given in terms of the effective coefficients (G- limits). 
In order to go further and design against failure initiation we record the effects of L°° constraints on 
the local gradient field inside functionally graded materials. For this we use the local representation 



formulas given by modulation functions (5.9) 



The multi-scale formulation of the graded material design problem has three features [23] , [24] : 
1. It admits a convenient local paramctrization of microstructural information expressed in terms 



ulation functions (5.9) 



of a homogenized coefficient matrix (5.6) and local representation formulas given by the mod- 



2. Is well posed, i.e., an optimal design exists. 

3. The optimal design is used to identify an explicit "functionally graded microstructure" that 
delivers an acceptable level of structural performance while controlling the local gradient field 
over a predetermined part of the structural domain. 

In what follows we describe the multiscale material design problem and focus the discussion on the 
control of the L°° norm of the local gradient field. 
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The admissible set of continuously graded locally periodic microstructures is specified by a vector 
P = (3 n ) of local geometric parameters. For example one may consider a periodic array 

of spheroids described by the orientation of their principle axis and aspect ratio. The periodic 
microstructure is specified in a unit period cell Y centered at the origin. Points in the cell are 
denoted by y. The characteristic function of the i th phase in the unit cell is denoted by X^/?,?/), 
i = 1, . . . , N. The vector f3 for the graded microgeometry can change across the structural domain 
SI and we write (3 — f3(x) for x 6 ft. The x dependence of f3 corresponds to the gradation of material 
properties through a gradation in microstructure. The design vector fi(x) is a uniformly Holder 
continuous function of x in the closure of f2. We write 

X l (x,y)= X l W(x),y) (6.1) 



and since fi(x) is continuous one sees that X l { x iV) i s continuous in the sense of (5.1). 

The multi-scale design problem is formulated as follows: The admissible set Ad of design vectors 
(3(x) is the set of uniformly Holder continuous functions satisfying the two conditions: 



There is a fixed positive constant C such that: 



\/3(x) -P(x')\ , s 

sup j-y^-; <c. (6.2) 

• The design vector (3{x) takes values inside the closed bounded set given by the constraints 

k < Pi{x) < b it i = l,...,n. (6.3) 

The local volume fraction of the i th phase in the composite is given by 8i(x) = J„ x l { x i v) dy- A 
resource constraint is placed on the amount of each phase appearing the design. It is given by 



6i(x)dx < 7i , i = l,...,N. (6.4) 
The vector of constraints (71, . . . , 7jv) is denoted by 7. The set of controls (3(x) £ Ad that satisfy 



the resource constraints (6.4) is denoted by Ad 1 . 

As an example we assume homogeneous Dirichlet conditions on the boundary of the design 
domain f2. For a given right hand side / € iJ -1 (f2) the overall structural performance of the graded 
composite is modeled using the solution u H of the homogenized equilibrium equation given by the 
i/g(£l) solution of 

-div (A H {x)Vu H ) = f. (6.5) 



Here A H is given by (5.6) with Xi{ x iV) given by (6.1) 



In this example the overall work done against the load is used as the performance measure of the 
graded material structure. This functional depends nonlinearly on the design j3 through the solution 
u H and is given by 

W(§) = [ fu H dx, (6.6) 
Jn 

We pick an open subset S C f2 of interest and the gradient constraint for the multi-scale problem 
is written in terms of the modulation function. We set 

C i (§) = \\M i (Vu H )\\ L o* is) , iori = l,...,N (6.7) 

and the multi-scale optimal design problem is given by 

P= inf {W(/3) : CM) < M, i = 1,...,N}. (6.8) 

peAd^ — — 
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When the constraint M is chosen such that there exists a control /3 € Ad 1 for which Ci{j3) < M 



then an optimal design /3 exists for the design problem (6.8), this is established in [28], [22 



The optimal design f3* specifies characteristic functions x l (x,y) = X l {P*( x )iy) from which we 
recover continuously graded microgeometries x%- (x,Tijx) and coefficient matrices 
A*' kj (x, rijx) of the form (5.5). The coefficients A* ,kj { x, rij x) G-converge to the effective coefficient 



A associated with the optimal design (3 , see Lemma 5.2 Here the effective coefficient is given by 



(5.6 ) with x 2 ( x j v) = X 1 { x i y)- F° r each j = 1, . . . the Hq (Q) solution Uj of the equilibrium problem 



inside the graded composite satisfies 

- div (A*' kj (x, njx)Vuj) = f 



(6.9) 



and the work done against the load is given by W(uj) = fujdx. This functional is continuous 
with respect to G-convergence hence lhrij_j. 00 W(uj) — W((3*). 

We now apply Theorem |2.2| to discover that for any open set S C £1 with closure contained inside 
there exists a decreasing sequence of sets for which \E^ j \ \ and 



limsup||xL(a;,njx)Vu i (x)||£oo (S \ B ) < M, i = 1,2, 



(6.10) 



Therefore we can choose a graded material design specified by Xk ( x J n j x ) with overall structural 
properties W(uj) close to the optimal one W(f3*) and with 



\\Xk j (x,n j x)Vu j (x)\\ L0 o {S \ Ek . ) 



< M 



outside controllably small sets E^- ■ This is the essence of the design scheme for continuously graded 
composite structures developed in [55], [2"^] . 

We conclude this section with a conjecture. Numerical simulations [24] show that when the 
microstructure corresponds to smooth inclusions embedded inside a matrix, such as shafts reinforced 
with long prismatic fibers with circular cross section, then the design method implies full control of 
the local gradient over the set S i.e., 



\\x\ A x ,n j x)'Vu j {x)\\ L ^ {s) < M, i = 1,2, 



,M. 



(6.11) 



With this in mind and in view of Theorem |4.2| we are motivated to propose the following conjecture. 

Conjecture 6.1. For continuously graded composites containing inclusions with C ' a boundaries 
for which A k i(x,nj) G-converges to A H (x) then 



limsupllxfc/^^i^VujIlioofg) = \\M l (Vu H )\\ L ^ (S )- 



(6.12) 



Appendix 



Here we will show that the solutions u>* of the cell problem (5.8) satisfying J Y w l (x,y)dy = are in 



C(f2, Hp er (Y)) under the continuity assumption (5.1 ). To that end, it suffices to show that as h — > 



\\V v w l {x + h,-)- V y w l {x, -)llL 2 (y) 0. 



Since w l (x + h,y) solves equation (5.8) A(x,y) replaced by A(x + h,y), we have that 



div (A(x + h, y){V v w l {x + h,y) + e 1 )) = div (A(x, y)(V y w t (x, y) + e 1 )) = 0. 
Rewriting the above equation we obtain 

div [(A(x + h,y)- A(x, y))}(\7 y w\x + h,y) + e l ) = div (A(x, y)(\7 y w l (x, y) - V y w l (x + h, y)). 
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Define the difference mapping Sh{F) = F(x + h,y) — F(x, y). Then for any ip £ Hp er (Y), we have 
(A(x,y)V y 5 h (w i )(x,y),ViP)dx= f (Sh(A)(x, y){V y w l (x + h,y) + e l ), Vip). (6.1) 



Then plugging if)(x,y) — Sh(w l )(x,y) £ Hp er (Y) in (6.1) and using the uniform ellipticity of the 
coefficients, we have 

X\\6 h (w l )(x, -)||i 2(y) < J (S h (A)(x, y)[V y w\x + h,y) + e% W^w^x, y))dy 

\ 1/2 

\S h (A)(x, y)[V y w\x + h,y) + e l ]\ 2 \ \\S h (w)(x, OlU^r) 



< 



The last inequality implies that 

AT , \ 1/2 

\\S h (w*)(x, OlU^y) < A/A^ / \ X l (x + h,y)- X \x, y)\ 2 \V y w\x +h,y) + e*\ 2 dy 

i=i ^ Y 

By Meyer's higher regularity result, VyW l (x + h, ■) £ L P (Y) for some p > 2. Moreover, the L p norm 
is bounded from above by a constant C independent of x, and h. After applying Holder's inequality 
we get 

\\6 h W)(x,-)\\ LHY) < — £ / |x i (» + fc,»)-x < (a;,»)| 2 di/ 

i=l ^ Y 



Applying (5.1 1, the right hand side approaches as h —> and the proof is complete. 



Acknowledgment 

This work is supported by grants: NSF DMS-0807265 and AFOSR FA9550-05-0008. Tadele Menge- 
sha gratefully acknowledges the support of Coastal Carolina University, where he is an assistant 
professor and currently on leave of absence. 

References 

[1] G. Allaire, Homogenization and two-scale convergence, SI AM J. Math. Anal, 23, (1992) 1482- 
1518. 

[2] M. Avellaneda, and F.-H. Lin, Compactness methods in the theory of homogenization, Comm. 
Pure Appl. Math. 40, no.6, (1987 ) 803-847 

[3] A. Bensoussan, J.L. Lions, and Papanicolaou, G. Asymptotic Analysis for Periodic Structures. 
Studies in Mathematics and its Applications, (5, North-Holland, Amsterdam, 1978 ) 

[4] B. V. Boyarskii, Generalized solutions of a system of differential equations of the first order of 
elliptic type with discontinuous coefficients, Mat. Sbornik N. S. 43, (1957) 451-503. 

[5] E. Bonnetier and M. Vogelius, An elliptic regularity result for a composite medium with "touch- 
ing" fibers of circular cross-section, SI AM J. Math. Anal. 31, (2000) 651-677. 

[6] L. A. Caffarclli and I. Peral, On W 1,p estimates for elliptic equations in divergence form, Comm. 
Pure Appl. Math. 51, (1998) pp. 1-21. 



24 



[7] J. Casado-Diaz, J. Couce-Calvo, and J. D. Martin-Gomez, "Relaxation of a control problem in 
the coefficients with a functional of quadratic growth in the gradient, SIAM J. Control Optim., 
47, ( 2008) 1428-1459. 

[8] M. Chipot, D. Kinderlehrer, and L. Vergara-Caffarelli, Smoothness of Linear Laminates. Arch. 
Rational Mech. Anal, 96, no 1 (1985) 81-96 

[9] E. De Giorgi, and S. Spagnolo, Sulla convergenza degli integrali dell' energia peroperatori ellittici 
del secondo ordine, Boll. UMI 8, (1973) 391-411. 

[10] P. Duysinx, and M. P. Bendsoc, Topology optimization of continuum structures with local stress 
constraints, Int. J Numer. Math. Engng., 43, (1998) 1453-1478. 

[11] J. Carlos-Bellido, A. Donoso, and P. Pedregal, Optimal design in conductivity under locally 
constrained heat flux. Arch. Rational Mech. Anal, 195, (2010) 333-351. 

[12] D. Faraco, Milton's conjecture on the regularity of solutions to isotropic equations Annales de 
L 'Institute Henri Poincare (c) Nonlinear Analysis 20, (2003) 889-909. 

[13] D. Fujii, B.C. Chen, and N. Kikuchi, Composite material design of two-dimensional structures 
using the homogenization design method, Internat. J. Numer. Methods Engrg., 50, (2001) 2031- 
2051. 

[14] D. Gilbarg and N.S. Trudinger, Elliptic partial differential equations of second order. (Springcr- 
Verlag, Berlin, New York, 2001). 

[15] S. Jimenez, and R. Lipton, Correctors and field fluctuations for the p e (a;)-Laplacian with rough 
exponents." J. Math. Anal. Appl. 372, (2010) 448-469. 

[16] J.H. Gosse, and S. Christensen, Strain invariant failure criteria for polymers in composite ma- 
terials, AIAA (2001) 2001-1184. 

[17] V. V. Jikov, S. M. Kozlov, and O. A. Oleinik, Homogenization of Differential Operators and 
Integral Functionals, (Springer- Verlag, Berlin, New York, 1994). 

[18] A. Kelly and N.H. Macmillan, Strong Solids. Monographs on the Physics and Chemistry of 
Materials. (Clarendon Press, Oxford, 1986). 

[19] Y.Y. Li, and M. Vogelius, Gradient estimates for solutions to divergence form elliptic equations 
with discontinuous coefficients, Arch. Rational Mech. Anal. 153, (2000) 91-151. 

[20] Y.Y. Li , and L. Nirenberg, Estimates for elliptic systems from composite material, Comm on 
Pure and Appl. Math. LVI, (2003) 892-925. 

[21] F. Leonetti, and V. Nesi, Quasiconformal solutions to certain first order systems and the proof 
of a conjecture of G.W. Milton, J. Math. Pures. Appl. 76, (1997) 109-124. 

[22] R. Lipton, and M. Stuebner, Inverse homogenization and design of microstructure for point wise 
stress control, Quarterly Journal of Mechanics and Applied Mathematics, 59, (2006) 139-161. 

[23] R. Lipton, and M. Stuebner, Optimal design of composite structures for strength and stiffness: 
an inverse homogenization approach, Structural and Multidisciplinary Optimization 33, (2007) 
351-362. 

[24] R. Lipton, and M. Stuebner, A new method for design of composite structures for strength and 
stiffness. 12th AIAA/ISSMO Multidisciplinary Analysis & Optimization Conference, Septem- 
ber 10-12, 2008, Victoria British Columbia, Canada. American Institute of Aeronautics and 
Astronautics Paper AIAA 2008-5986. 



25 



[25] R. Lipton, Assessment of the local stress state through macroscopic variables, R. Soc. Lond. 
Philos. Trans. Ser. A Math. Phys. Eng. Sci., 361, (2003) 921 -946. 

[26] R. Lipton , Bounds on the distribution of extreme values for the stress in composite materials, 
J. Mech. Phys. Solids, 52, (2004) 1053-1069. 

[27] R. Lipton, Homogenization and field concentrations in heterogeneous media. SIAM J. on Math- 
ematical Analysis, 38, (2006) 1048-1059. 

[28] R. Lipton, Homogenization and design of functionally graded composites for stiffness and 
strength, In P. Ponte Castaneda et. al. eds., Nonlinear Homogenization and its Applications 
to Composites, Polycrystals and Smart Materials, Kluwer Academic Publishers, Netherlands, 
(2004), 169-192. 

[29] A.J. Markworth, K.S. Ramesh, and W.P. Parks, Modelling studies applied to functionally graded 
materials, Journal of Materials Science 30, (1995) 2183-2193. 

[30] N. Meyers, An L p -Estimate for the gradient of solutions of second order elliptic divergence 
equations, Annali della Scuola Norm. Sup. Pisa 17 (1963), 189-206. 

[31] G.W. Milton, Modeling the properties of composites by laminates, In: Homogenization and 
Effective Moduli of Materials and Media Edited by J. Erickson, D. Kinderlchcr, R.V. Kohn, 
and J.L. Lions. IMA Volumes in Mathematics and Its Applications 1, (1986) 150-174. Springer- 
Verlag, New York. 

[32] F. Murat, and L.Tartar, Calcul des Variations et Homogeneisation, Les Methodes de 
I'Homogeneisation: Theorie et Applications en Physique, Edited by D. Bergman et. al., Col- 
lection de la Direction des Etudes et Recherches d'Electricite de France, 57, Eyrolles, Paris, 
(1985) 319-369. 

[33] G. Nguetseng, A general convergence result for a functional related to the theory of homoge- 
nization, SIAM J. Math. Anal., 20, (1989) 608-623. 

[34] R.J. Nuismer and J.M. Whitney, Uniaxial failure of composite laminates containing stress con- 
centrations, In: Fracture Mechanics of Composites. ASTM Special Technical Publication Amer- 
ican Society for Testing and Materials, Philadelphia. 593, (1975) 117 142. 

[35] Y. Ootao, Y. Tanigawa, and O. Ishimaru, Optimization of material composition of function- 
ally graded plate for thermal stress relaxation using a genetic algorithim, Journal of Thermal 
Stresses, 23, (2000) 257-271. 

[36] G. Papanicolaou, and S.R.S. Varadhan, Boundary value problems with rapidly oscillating ran- 
dom coefficients, Random fields. Rigorous results in statistical mechanics and quantum field 
theory, Esztergom 1979, Colloq. Math. Soc. Janos Bolyai, 27, (1981), 835-873. 

[37] E. Sanchez-Palencia, Non Homogeneous Media and Vibration Theory (Springer, Heidelberg, 
1980). 

[38] S. Spagnolo, Convergence in Energy for Elliptic Operators. Proceedings of the Third Symposium 
on Numerical Solutions of Partial Differential Equations, Edited by B. Hubbard, (College Park, 
1975). Academic Press, New York, (1976) 469-498. 



26 



